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Abstract 

The FAA is currently developing the Wide Area Augmentation System (WAAS), a GPS Satellite-Based 
Augmentation System (SBAS) that will be used for navigation within the National Airspace System 
(NAS). WAAS consists of a multitude of ground reference stations and several geo-stationary 
communication satellites. The geo-stationary satellites will be broadcasting a “GPS”-like signal that can 
be used as a supplementary navigational signal and that will also contain GPS health status and 
correction data. The WAAS also offers a significant tool to the timing community. 

There are several unique advantages of using the WAAS or any SBASfor frequency and time transfer: 

1. The geo-stationary satellites are always in view. They never set like the GPS and GLONASS 
satellites. This offers the ability for multiple sights to be permanently “phase-locked” together 
by locking themselves to the common WAAS signal. 

2. They do not have any intentional distortions like the GPS SA dither. 

3. WAAS system computes and broadcasts orbit and clock corrections for all GPS satellites in view 
thereby substantially removing the effects of SA on GPS. 

4. The geo-broadcast signals are cesium-based, generated and controlled by the WAAS network of 
cesium-based GPS receivers. 

5. The WAAS network computes a real-time observational ionospheric delay model for the 
coverage area. This model is broadcast over the WAAS signal so that an LI single frequency 
user can make local corrections. 

This paper describes an initial attempt at time transfer using the WAAS test transmissions currently 
broadcast from the INMARSAT Atlantic Ocean Region - West (AOR-W) satellite. The clocks being 
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compared are the master clock at USNO Washington and the master clock at NRC Ottawa. A 30-day 
data set is examined. It was obtained from dual frequency GPS/WAAS NovAtel receivers that were set 
up at both institutions. Comparisons are made between UTC(USNO), UTC(NRC), GPS time,and WAAS 
Network Time (WNT). 

I. Introduction 

A NovAtel WAAS MiLLennium Narrow Band Receiver was installed and operated at the United 
States Naval Observatory's (USNO). The receiver's oscillator was phase-locked to the USNO 
hydrogen maser master clock (MC2). The location of the temporary antenna installed on the roof 
was computed from 3 weeks of averaging of WAAS corrected GPS computed points. There were 
107,420 points gathered. One point was collected every 15 seconds. The average computed 
position was: 

N 38° 55’ 13.8138", 

W 77° 03’ 59.2032", 

Height 53.044 m (Ellipsoid) 

The standard deviation of each point. Latitude, Longitude, and Height was 1.3 m, 1.2 m and 3.2 m 
respectively. 

A second NovAtel MiLLennium receiver was installed at the Canadian National Research 
Council (NRC) in Ottawa. It was remotely operated from Calgary over the Internet with 
RapidRemote VI.5 software from QuarterDeck. This receiver was connected and phase-locked 
to NRC's hydrogen maser master clock (H4). The coordinates of NRC were also computed by 
taking a 3-week average. The average computed position was: 

N 45° 27'14.9814" 

W 75° 37' 25.7839" 

H 82.732m (Ell.) 

The precise vector between NRC and USNO was computed by processing three separate 24-hour 
data sets using NovAtel's SoftServe Geodetic dual frequency postprocessing package. Holding 
NRC fixed at its WAAS averaged computed location, the phase vector position computed for 
USNO was: 

N38° 55' 13.7959" 

W77° 03" 59.1977" 

H 51.784m (Ell.) 

The vector distance between the two receivers was 759Km on a north-north-east baseline from 
USNO. The agreement between the three successive processing batches was better than 5cm. The 
difference between the 3-week averaged USNO point and the computed vector from NRC 
was 0.5m, 0.1m, 1.3m for Latitude, Longitude, and Height respectively. The vector-computed 
location for USNO was used for all processing because it was important to maintain a very 
accurate relative position between USNO and NRC for the very clock transfer computations. 
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II. Data Collection 


Data collection was commenced Oct 6 1999 and continued through to Nov 31 1999. Pseudo 
range and carrier phase measurements were made at 15-second intervals from GPS as well as the 
AOR-W'WAAS GEO satellite. All of the broadcast data from the satellites was collected for 
subsequent postprocessing. Most of the analysis for this paper used 3 weeks of data from >m 
November 6 through to November 22. There were several data breaks and setup problems with 
the October data. 


III. Data Processing 


The clock difference between any sight (say USNO) and GPS system time was computed by 
Equation (1): 


UPSLlockdff] (/) = Psr(t) i + Trope ft) + lonocft ) - RangeloSatft) (1) 

Where: 
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is the time of measurement 

is the clock difference at time t between the receiver’s time and the GPS 
system time as measured by satellite i. 
is the measured pseudorange to satellite i 
is the Tropospheric correction for satellite / 
is the Ionospheric correction for satellite / 

is the computed range to satellite i using GPS broadcast ephemeris 


GPSClockdiff (t) = —^—y GPSCiockdifft) (2) 
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Note that the receiver clock difference at any point in time can be further improved by averaging 
the computed clock differences from ail satellites observed at that time using Equation (2). 
Typically there are around 8 satellites in view with a minimum of 4 and a maximum of 12. The 
NovAtel Millennium WAAS receiver was configured as 10 GPS L1+L2 plus 1 WAAS LI. This 
allowed for a maximum of 11 satellites tracked at any point in time. 


To analyze the frequency stability of the Clockdiff measurement, the Allan deviation curve was 
computed using Equation (3). 


AllanDeviation(T) = 
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Where: 

T is sample period (s) 

N is number of sample sets used in accumulation 


To analyze the time stability of the Clockdiff measurement, a 2nd order low pass feed back filter 
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was used to "smooth" the data. The design of this "Smoother" filter was chosen to closely match 
a circuit that would slave a local GPS oscillator that was producing the Clockdiff measurement. 
Fig. 1 illustrates this filter. 

The Clockdiff measurements are processed by passing the data repeatedly through the smoother 
filter with an adjustable time constant. The standard deviation of the smoothed data is computed 
for each time constant used. The low pass smoother characteristics indicate whether an 
improvement in timing accuracy can be made by averaging the instantaneous Clockdiff 
measurements with the use of a good local oscillator. 


Frequency transfer function for Low pass Smoother second order loop is given by Equation (4): 
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The clock difference between any receiver and WAAS Network Time (WNT) was 
computed by Equations (11) and (12): 

WNTClockdiff, = Psr, + Trope , + Ionoc, + WAASFast, + WAASSlow>, - RangeToSat , ( 1! ) 
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Where: 


WNTClockdiffi is the clock difference between the receiver's time and the WNT as 
measured by satellite /. 

WAASFasti is the fast clock correction for satellite i 

WAASSloWi is the slow orbital error and slow clock correction for satellite i 

The WNTClockdiff was analyzed by computing its Allan deviation and smoother characteristics 
by Equations (3) and (4). 

The clock difference between two sights on the ground (say USNO and NRC) can be computed 
using either Equation (13) or (14). 

DClockdiff, = GPSClockdiff™ - GPSClockdiff, mc (l3) 

DClockdiff, = WNTClockdiff, USNO - WNTClockdiff (14) 

Where DClockdiff is the clock difference between the USNO and NRC clocks as measured by 
satellite i. 

It would be expected that the Dclockdiff computed from Equation (14) using the improved 
WAAS orbit data would significantly improve the results as compared with only the broadcast 
GPS positions as computed by Equation (13). Equation (14) was used for the results of this 
paper. 

The correction for the ionosphere (Ionoci) was computed from the difference between GPS 
smoothed LI and L2 pseudorange measurements.Fig. 2 illustrates the ionospheric correction 
filtering for a single GPS L1-L2 satellite pass. The differenced pseudorange is smoothed by the 
differenced carrier phase measurement. Essentially, the differenced carrier phase values are used 
for the ionosphere after they are adjusted to the average offset that maps it onto the differenced 
pseudorange curve. 

The pseudorange value, corrected for ionospheric delay, was similarly represented by the LI 
Carrier measurement, corrected for ionospheric delay. The carrier provided a low noise estimate 
of the pseudorange value. Fig. 3 illustrates this process. Fig. 3 shows the misclosure between the 
carrier phase and the pseudorange value. Note that the carrier phase fits the pseudorange 
measurements with zero mean and 0.1 ns standard deviation over this 6-hour pass. Most of the 
excursions between the code and carrier are caused by multipath interference on the pseudorange 
measurement at the receiving antenna. 

Only GPS measurements from satellite passes that were longer than 3 hours where the satellite 
was 20 degrees above the horizon were used in the calculations. Having a high elevation cut-off 
angle removed most of the effects of multipath interference in the data. Having long sets of 
satellite data enabled very accurate carrier fit estimates for the ionospheric correction and the LI 
pseudorange smoother. 

A standard Hopfield model was used for the tropospheric correction. 
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IV. Discussion of Results 


The difference between USNO and GPS system time was computed and is illustrated in Fig. 4. 
This difference was included as a common denominator benchmark process. Most GPS timing 
receivers today use this observable to provide time reference information. Note that the 43ns 
standard deviation caused by the DOD's selective availability (SA) (intentional dither) is very 
slow in nature and must be smoothed over a very long time constant in order to be removed. The 
"Smoother" curve in Fig. 5 illustrates the correlation characteristic of this process. Note, in Fig. 5 
the smoothing time constant must be at least 55,000 seconds before a significant filtering of the 
SA is accomplished and has been reduce to less than 10ns. Note also that the Allan deviation is 
dominated by SA starting out at 4x1 O' 10 at a 30-second measurement interval and falls at a 
constant lOdB per decade, showing that the frequency dither of SA is uncorrelated. 

Fig. 6 represents a 3-week history of the difference between the clock at USNO and WNT. The 
standard deviation of the clock difference during this period was 11 ns. But there was significant 
low frequency noise presumably caused by clock steering oscillations in the WAAS system 
control segment. We make this assumption because the oscillations occur on even 12-hour 
boundaries (coincident with the WAAS steering algorithm). Fig. 7 shows the difference between 
the USNO clock and WNT during the later 12 hours of November 14th. Note that the fit between 
USNO and WNT over short 12-hour segments is in the nanosecond noise level. This result is 
very significant in that once the WAAS control segment smoothes out its clock steering 
oscillations, the long-term clock errors will be at the nanosecond level. 

Fig. 8 shows the Allan deviation and smoothing characteristics of the 3-week USNO-WNT data 
set. Note, in Fig. 8, that the Allan deviation flattens out to 3x1 O' 13 at 10,000 seconds because of 
the network steering oscillations. Similarly, the smoother time constant does not significantly 
filter the oscillations until its time constant is over 2 days in length. 

Fig. 9 shows the differential time drift between USNO and the NRC master clock. We used the 
WAAS corrected orbit Equation (14) to produce these results. The standard deviation from a 
linear fit (of 8ns/day) was less than Ins! 

Using the GPS broadcast orbits increase the standard deviation of the clock difference by a factor 
of 2, (an increase from 0.6 to 1.2 ns). 

Fig. 10 shows the Allan deviation and smoothing characteristics of the WAAS orbit differential 
timing data. It shows that the smoother does not converge until the time constant gets longer than 
the longest satellite pass tracking time (~10,000s). This is caused by the correlation of the phase 
smoothing process of the ionosphere and the LI Pseudorange. Note that the standard deviation of 
the Smoothed clock offset drops below 400ps after 1 day of smoothing. 

The Allan variance of the differential data starts off at 6x1 O' 13 at 120-second sample period and 
reduces linearly down to 3E-15 at 2.5-Day sample period (200,000s). 
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V. Conclusions 


The receiver was able to track WNT with a noise level of 1.5ns. Using the WAAS computed 
orbits provided the common view multiple satellite time transfer of less than Ins (measured 
0.7ns). This method of providing time transfer is simple and effective. The biggest sources of 
error were the transitions between satellites (constellation switches). 

These results, at this early stage in the development of the WAAS, indicate the future promise of 
this technique for time transfer and time distribution. 
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Questions and Answers 


HUGO FRUEHAUF (Odetics): How does the fact that the GUS is, of course, stationary and 
the satellite is sort of making a figure-eight ground trace — is there a correlation in the timing 
data that you see from that figure eight? 

PATRICK-FENTON (NovAtel): Well, the figure eight is plotting a very slow multi-path that 
we see in the diagram. The multi-path, instead of in a typical GPS multi-path, is in the order 
of 3 to 5 minutes, where the GUS multi-path is like 5 to 6 hours. So you can see that errors 
correlated to the little figure eight. But as far as using the geostationary satellites, it’s really 
nice because it’s not moving and you can put a dish on it and get around that kind of stuff. 
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